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Abstract 

The nonlinear propagation of low-frequency circularly polarized waves in a magnetized dusty 
plasma is analyzed. It is found that wave steepening and shock formation can take place due to 
the presence of nonlinear quantum vacuum effects, thus giving rise to ultra-intense electromagnetic 
shocks. Moreover, it is shown that solitary wave structures are admitted even under moderate 
astrophysical conditions. The results may have applications to astrophysical plasmas, as well as 
next generation laser interactions with laboratory plasmas containing dust clusters. 

PACS numbers: 52.27.Fp, 52.35.Mw, 52.38.-r, 52.40.Db 
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I. INTRODUCTION 



The extreme environment offered in many astrophysical situations, such as supernovae 
and binary mergers, in many cases give rise to new and exciting physics, e.g. neutrino-plasma 
interactions and gravitational waves. Furthermore, the next generation laser-plasma systems 
promises to take laboratory astrophysics into new realms, as tens to hundreds of petawatts, 
or even zetawatts, may be generated in femtosecond pulses Q|. Thus, as observational 
astronomy and cosmology grows even more accurate, and laboratory techniques continue 
to develop at the current rate, one will be able to probe new physics, such as neutron star 
equations of state Q or the Unruh-Hawking effect 

Recently, it has been noted that the inclusion of quantum electrodynamics (QED) into 
Dlasma theory may give rise to novel wave modes J, S, 1^ , in particular in dusty plasmas 
7|. These wave modes could propagate in a variety of environments, such as supernovae 
remnants, interstellar clouds, or laser-dusty plasma systems. 

In Ref. jsl it was found, using the method of characteristics, that ultra-strong magnetic 
fields may induce shock front formation in electromagnetic radiation, analogous to hydrody- 
namic shocks. A possible application to such a scenario is magnetars, where field strengths 
of the order 10^^ G can be found. 

In the present paper, we show that the nonlinear propagation of electromagnetic waves in 
a a weakly magnetized dusty plasma may give rise to wave steepening and shock formation. 
This effect is solely due to the novel QED effect of photon-photon scattering in a nonlinear 
quantum vacuum. The shock front formation may give rise to ponderomotive acceleration 
of charged particles in, e.g. interstellar clouds. The results are discussed for astrophysical 
energy scales. 



II. PULSE EVOLUTION 



In Ref. m the dispersion relation 
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was found. The relation (^J) is valid for low-frequency circularly electromagnetic polarized 
waves, with electric field E = £'(x + iy) ex\i{ikz — iut), propagating along a magnetic 
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field Bq = Bqz Moreover, the critical Schwinger field is given by Es = m^^c^/eTi ~ 10^^ 



V/m (see, e.g. Ref. where c is the speed of light in vacuum, Ti is Planck's constant 
divided by 27r, moe is the electron rest mass, e is the absolute value of the electron charge, 
and ujcd = qdBo/mQd and cOpd = (^od^J/^o^^od)"'^^^ is the dust gyro and plasma frequency, 
respectively. Here uq^ denotes the dust particle density in the laboratory frame, the dust 
charge, mod the dust particle rest mass, and a is the fine-structure constant. 

Next, we introduce the dimensionless variables and parameters Q = uj/ujcd-, K = kc/ujpd, 
and K = {Aa / AbnY^'^ujpd/ ujcd, and normalize the electric field according to £" = E/Eg. We 
consider the case of ultra-relativistic dust in a weak magnetic field, for which 0| 

K^n^ + n^ = K'\£\^K\ (2) 

holds. 

The generic evolution equation for a slowly varying envelope £ following from the disper- 
sion relation (jSJ reads 

idl8-dldl£ + K^\£\^dt£ = Q, (3) 

where iK dz and iVt —dt has been used. 

If the quantum parameter k — 0, from Eq. (j2^ve obtain the equation 

{idt-dl)£(t,z) = Q, (4) 

yielding a linear traveling wave solution with the whistler-like dispersion relation Vt = —K^ 
In the case of small frequencies, Vt ^ K^, (the case with a dominant nonlinear term), we 
neglect the third order frequency VL^ in Eq. (2) and obtain 

dt£±K\£\dz£ = (5) 

Thus, letting £(t,z) = \£{t, z)\exp[iip{t, z)] the amplitude |£^(t,2)| satisfies the equation 

dt\£\±K\£\dz\£\ = 0, (6) 

while the phase is connected with the amplitude by the equation dt^p ± K\£\dz(p = 0. For 
simple waves 13] Eq. (jH)) yields the solution 

z = ±K\£\ + f{\£\), (7a) 
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I.e. 

\£it,z)\=f-\zTf^m, (7b) 

where the function is the inverse function of / and / is determined by the initial condition 
at t = 0. Thus, we have to define the initial condition |£^(0,z)| = f{z) at t = and then 
find z = f^^{\S\). As Eq. ()7b|) shows, the wave steepening due to nonlinearity may occur. 
This can also be seen via Eq. ()7a|). from which we find that 



where f'{x) = dxf{x) 



(a) If f'{\S\) < (which can be interpreted as the initial function \S\ = f{z) at t = 
decreasing with z), we have at time t = —\f'{\S\)\/K that dz\S\ = oo, if the plus is 
chosen in front of the K-term in Eq. (jHJ. 

(b) If we, on the other hand, have f'{\S\) > 0, we choose the minus sign in front of the 
K-term in Eq. (jHl), and obtain the shock front formation for a finite time. 

Thus, we clearly have the shock wave formation at this instant .lOj . Starting with a Gaussian 
pulse, case (a) implies the wave steepening and the shock front formation in the leading part 
of the pulse, while case (b) implies the wave steepening in the tail of the pulse, which can 
be seen by a simple spatial reflection symmetry. We note that in principle density gradient 
effects in the dusty plasma should be taken into account once significant steepening occurs, 
but this turns out to be a higher order effect. 

In the weakly nonlinear case, we may use Eq. (@]) to simplify the general pulse evolution 
equation 0, using d'^S ~ idtS. Assuming a slowly varying envelope S, Eq. yields 

idtS - djS - K^\S\^S = 0, (9) 

which is the well-known nonlinear Schrodinger equation with a focusing nonlinearity. 
Thus, such an equation admits the formation of a bright soliton given by S{t,z) = 
£^0 sech(£^o-2) exp (— z£^Qt/2v^) , where S'q is the amplitude, and we have assumed that the 
soliton has the peak value at t = located at 2; = 0. Thus, in the weakly nonlinear case, a 
balance between the dispersive and nonlinear effects may be struck. However, in practice, 
the coherence length of such structures could be short, due to e.g. plasma inhomogeneities. 
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Moreover, while the one-dimensional bright soliton is stable, higher dimensional effects may 
yield instabilities. It is interesting thought that such structures could form with rather 
modest parameter values, as will be discussed next. 



III. DISCUSSION AND CONCLUSIONS 

The effect of the wave steepening and intense shock front formation presented here takes 
place whenever the electromagnetic intensity is high in the dusty plasma. Thus, a natural 
question arises: What are the situations when this may occur? As astrophysical systems 
offer extreme energy scales, these are natural places to look for QED effects, such as the one 



le presence of dusty plasmas is we 



and interstellar clouds 



1 known in a wide range of 



presented here. Moreover, t 
astrophysical environments 

In the case of Interstellar clouds, we have an expected dust density of the order 
nod ~ 10~^ m~^ [17], with a typical dust mass mod ~ lO^mproton ~ 10"^^ kg [l^. A value 
Bq ~ 10^*^ G (well within the limits for such environments) results in k ~ 10^. If the 
circularly polarized wave satisfy Eq/Es ~ 2 x 10~^, a value appropriate at a distance of 
lAU from a supernova, we have k\S\ ~ 1, thus giving an appreciable nonlinearity in Eqs. 
(0) and ©. 

In Fig. 1 we have plotted the pulse evolution in two different cases. The upper panel 
represents the positive sign of Eq. (jH)), while the lower panel depicts the pulse evolution 
for the negative sign in Eq. ®. We have used the initial pulse amplitude |£^(0,2;)| = 
£oexp[—{z — zq)"^ /2o?], with zq = ?> and a = 1. Moreover, the initial amplitude satisfies 
kEq = 0.7, thus well within the limits posed by, e.g. interstellar clouds at lAU from a 
supernova. We have plotted the evolution of \S{t, z)\/So. The effect of the nonlinearity can 
clearly be seen in the form of the wave steepening. 

One interesting aspect of the shock front formation is the possibility of charged particle 
acceleration due to the ponderomotive force. However, before conclusions can be drawn on 
possible particle acceleration, analysis of mechanisms rendering the pulse wave front less 
steep must be made. 
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FIG. 1: The magnitude \£\ of the normahzed electric field relative to the initial value £q plotted 
as a function of z for three different times. At t = we have the dashed Gaussian curve, while 
at i = 12 and t = 23 we have the dotted and full curve, respectively. In the upper panel we have 
used the positive sign in Eq. @. The wave steepening and shock front formation can clearly be 
seen. In the lower panel, we have instead used the negative sign in Eq. Q. As expected, the wave 
steepening and shock front formation now occurs at the trailing edge of the pulse. 
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